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Data obtained in space, particularly that from IMP-1, have estab- 
lished that many of the gross features of the interaction of the solar 
wind and the geomagnetic field can be described by the continuum theory 
of fluid flow. The purpose of this paper is twofold: (a) to indicate 
briefly how one can arrive at results for the magnetosphere boundary, 
the bow wave, and the associated flow field from a purely fluid point 
of view without resort to a mixture of particle and fluid models as 
usually employed and (b) to present numerical results for the flow param- 
eters in the shock layer between the bow wave and the magnetosphere 
boundary. The objective is not to get a complicated theory that encom- 
passes all phenomena of Importance, but rather the simplest body of anal- 
ysis that appears able to describe the average bulk properties of the 
solar wind as it flows steadily through the bow wave and around the for- 
ward portion of the magnetosphere. A more complete account of this study 
is given by Spreiter et. al . (l965)* 

The fundamental assumption is that the flow can be described ade- 
quately by the standard magnetohydrodynamic equations for steady flow of 
a dissipationless perfect gas. Wo attempt will be made here to justify 
this assumption in detail, except to remark that the presence of a weali 
aiid irregular magnetic field in the incident solar wind plasma seems 
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sufficient to couple the motions of the particles even in the absence of 
collisions. The differential equations are thus as follows (see^ e.g., 
Landau and Lifshltz, 196O) . 

V . pv 0 "I 


p(v . v)v + VP - - 1^ H X curl H = - ^ VH^ + (H . v)H 


Lit 


) ( 1 ) 




(y . y)S = 0 , p = 


where p, p, S, and v refer to the density, pressure, entropy, and veloc- 
ity of the gas, H refers to the magnetic field, 7 = Cp/cy, and Cp and 
Cy are constants representing the specific heats at constant pressure 
and constant volume. Important auxiliary relations for 7, temperature 
T, speed of sound a, internal energy e, and enthalpy h are 


7 = (W + 2)/N , p = pRT/p = nkT 
e cvT, h = cpT = e + p/p 


( 2 ) 


where W represents the number of degrees of freedom, R = (cp - cy)p 
= 8 . 31^x10”^ erg/°K, p - mean molecular weight = l/2 for fully ionized 
hydrogen plasma, n = numher of particles/cm^ = 2np where np is the 
nunher of protons/cm^, and k = Boltzmann’s constant - 1.38x10”^® erg/°K. 
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Although only first derivatives appear in this system of equations, 
the neglected dissipative terms are described by second derivatives. 

Their neglect requires that the gradients be small. In magnetohydrody- 
namics, as in gasdynamics, however, compressions tend to coalesce and 
steepen into shock waves of such small thickness that they can be consid- 
ered as virtual discontinuities. In addition, attraction between like 

currents tends to cause current distributions to collapse into thin 
sheaths, across which the magnetic field can be considered in the same 
sense to be nearly discontinuous. Mathematically, continuous solutions 

of the dissipationless differential equations cease to exist, and the 
flow is no longer governed solely by the equations given above. Mass, 
momentum, magnetic flux, and energy must still be conserved, however, 
and the following relations must hold between quantities on the two sides 
of any such discontinuity; 


[pv^] = 0 


1 


/ 




P^nV + i P + j n 


PMt 






^ 0 


- 0 , [Hn] - 0 


\ ( 3 ) 


P'^n (2 "^n^ '*■ ’'^n hjj " Hn^ • S 


= 0 


The subscripts n and t refer to components normal and tangential to the 
discontinuity surface and [Q] = Qi - Qq where subscripts 0 and 1 refer 
to conditions on the upstream and downstream sides of the discontinuity. 
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Five classes of discontinuities sire described by equation (3). Fol- 
lowing Landau and Lifshitz (1960)^ discontinuities which lie along stream 
lines (vn = 0) are called tangential discontinuities or contact discon- 
tinuities according to whether or not the normal component of the mag- 
netic field Hn vanishes. Discontinuities across which there is flow 
(vq / 0) are divided into three categories called rotational disconti- 
nuities^ and fast and slow shock waves. Some properties that distinguish 
the various discontinuities are summarized below: 

Tangential discontinuities 

Vn Hji =■ 0 ; [vt] / 0 , [gt] ^ ^ , [p] ^ 0 , [p + =■ 0 

Contact discontinuities 

Vn = 0 , Hn / 0 : [v] = [H] = [p] = 0 , [p] ^ 0 (5) 

Rotational discontinuities 

v^ = ±Hn/(Wp)^^'^ , [vtl = [Ht]/(4:xp)^^^, ' 

( 6 ) 

[p] = [p] ^ [vn] - [Hn] - [H^] 0 


Shock waves^ fast and slow 


vn / 0 , [p] > 0 , [Hn] = 0 , 


^P""n^fast ^ ^P^n)rot ^ Hn(p/4rt)^^^ > (pvn)gn 


slow 


(7) 


H-5 and H‘ 


2 / increase \ 
'^decrease J 


f fast 
yslow^^ 


through 


shock waves 
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Of i.liese, on].y the tangential discontinuity has properties compat- 
Ihlo with those required to describe -the boundary of the geomagnetic 
field. As in the classical Chapman-Ferraro theory based on particle 
concepts^ the condition = 0 holds and requires that there is no 
connectivity between the geomagnetic and interplanetary fields. Slight 
differences stem from the additional assumptions in the Chapman-Ferraro 
theory that the incident plasma is free of magnetic field and the outer 
magnetosphere is free of plasma so that j>o ^ Although neither 

of these statements is strictly true, estimates of the magnitudes of the 
gas pressure p aid the magnetic pressure lead to the conclusion 

that p « in the outer magnetosphere and that p » in the 

shock layer. So far as the flow outside the magnetosphere is concerned, 
the discontinuities at the magnetosphere boundary may thus be approxi- 
mated satisfactorily by those of the limiting case of a tangential dis- 
continuity in which there is a vacuum on one side, and negligible mag- 
netic pressure on the other side. It is demonstrated further by 
Spreiter et . al . ( 1965 ) that the gas pressure on the magnetosphere 
boundary is adequately approximated by the simple Newtonian formula 
p ^ Kp Vo^^cos^f' where the subscript 00 refers to the values in the 
undisturbed Incident stream, \|/ is the angle between the normal to the 
boundary and the velocity vector of the undisturbed stream, and K is 
a constant usually taken to be unity but which yields better accuracy if 
equated to 0.84 for 7-2. and 0.88 for 7 = 5/3- In this way, the cal- 
culation of the shape of the nmgneto sphere boundary is decoupled from 
the detailed analysis of the surrounding flow and reduced to the identi- 
cal mathematical problem described by the classical Chapman-Ferraro theory. 
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Examination of conditions typicil of the solar wind shows that the 
mass flux through the nose of ohe bow wave is characteristically 

about an order of magnitude greater than that through a rotational dis- 
continuity. The bow wave must therefore be a fast shock wave. 

Although it has been necessary to invoke the presence of a magnetic 
field to provide a mechanism for providing interactions between the par- 
ticles^ and to consider the complications of hydromagnetic theory to pro- 
vide an adequate description of the magnetosphere boundary as a tangen- 
tial discontinuity and the bow wave as a fast hydromagnetic shock wave, 
the magnetic field strengths typically encountered in the solar wind are 
sufficiently small as otherwise to produce only secondary effects on the 
flow. On this basis, one can justify the neglect of all the terms con- 
taining the magnetic field H in equations (l) and (3) when calculating 
the flow. The fluid motion then calculated is just exactly that indi- 
cated by ordinary gasdynamic theory. Deformations of the magnetic field 
are determined subsequently using the equations 

curl (H X v) - 0 , div H = 0 , [Hnl =0 (8) 

which are frequently interpreted as indicating the magnetic field lines 
move with the fluid. 

In this way, the exterior magnetosphere flow problem is reduced to 
a purely gasdynamic problem of supersonic flow past a given body, the 
shape of which is determined by solving the standard Chapman -Ferraro 
problem. Difficulties associated with the nonlinear and mixed elliptic- 
hyperbolic character of the governing partial differential equations 
are such, however, that it is necessary at the present time to approximate 
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'.ho maf'iieiof.pliere boundary with an axisyinmetric shape. The precise shape 
coLoct.ed Jo. ihat obtained by rotating the equatorial trace of the bound- 
ary given by the approximate solution of Spreiter and Briggs (1962) 
about an axis through the center of the earth that extends parallel to 
the velocity vector of the undisturbed incident stream. 

The principal gasdynamlc result presented by Spreiter and Jones 
(1963) was a plot of the calculated position of the bow shock wave asso- 
ciated with the simplified axisymmetric magnetosphere described in the 
preceding paragraph. The calculations were performed for a ratio of 
specific heats / of 2 ., and a free-stream Mach number of 8.7I. The 
latter was identified with the free-stream Alfv^n Mach number associ- 
ated with a representative choice of values for the density, velocity, 
and magnetic field of the incident solar wind. It is implicit in the 
analysis presented above, however, that the gasdynamic Mach number should 
more properly be identified with the free-stream Mach number = Woo/aco 
than the Alfv6n Mach number M^ = Vog(4jtp/H^)^^^ if jq is much greater 
than unity, as is indeed generally the case in the present applications. 
Since a value of 8.7I is also a reasonable choice for Mx>^ we turn now 
to a presentation of further details of the flow field for the same set 
of conditions as employed by Spreiter and Jones (1963). 

Figure 1 shows a plot of the magnetosphere boundary and shock wave 
position for 7=2 and = 8.7I in terms of a length scale in which 

the distance D from the center of the earth to the magnetosphere nose 
is unity. D is generally of the order of 10 earth radii, and fluctuates 


in response to variations in the incident stream in accordance with the 


expression ® - ^eHp_ V(atKp V where ag = 6.37x10® cm is the 


radiuK of the earth and Hp^ - 0.312 gauss is the intensity of the earth’s 
magnetic field at the equator. Also included on this figure are several 
solid lines representing streamlines, and broken lines representing 
characteristic or Mach lines of the flow. The latter correspond to 
standing compression or expansion waves of infinitesimal amplitude. 

They cross the streamlines at such angles that the local velocity com- 
ponent normal to the wave is exactly equal to the local sound speed. 

Mach lines thus exist only where the flow is supersonic; they are absent 
from the vicinity of the magnetosphere nose because the flow there is 
subsonic. 

Contour maps showing lines of constant density, velocity, temper- 
ature, and mass flux, each normalized hy dividing hy the corresponding 
quantity in the incident stream are presented in Figure 2. The results 
show that the density ratio remains close to the maximum value 

(7 + l )/(7 - 1 ) = 3 lor a strong shock wave in a gas with 7 - 2 . along 
nearly the entire length of the portion of the how wave shown. The gas 
undergoes a small additional compression as it approaches the stagnation 
point at the magnetosphere nose and then expands to less than free-stream 
density as it flows around the magnetosphere. The velocity remains less 
than in the free stream, however, throughout the same region. The temper- 
ature T/Tqo is closely related to the velocity ratio through the 
expression 




= 1 + 


(7 - 1)Moo'' 


1 - 




(9) 


derived from equations (l) and (2) with H - 0. Of particular interest 
is the tremendous increase in temperature of the solar wind as it passes 
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thnnin;]! the liow chock wave. If, for exang^le, the temperature of the 
inci-clcni, colar wind Ic 50,000° K, the temperature at the magnetosphere 
nose :3 s Indicated to he nearly 2,000,000° K. This value is consistent 
with the temperature possessed hy the gas in the solar corona before 
it is accelerated to the high velocities characteristic of the solar 
wind in the vicinity of the earth’s orbit. 

As noted above, calculation of the deformation of the magnetic 
field in the flow around the magnetosphere can be accomplished directly, 
once the flow field is determined, either by direct integration of 
eq^uation (8) or by considering the field lines to move with the fluid. 
Although the field lines are, in general, spatial curves, simplicity 
may be achieved at the expense of completeness by confining attention to 
the plane containing the velocity and magnetic field vectors in the 
Incident stream. Since the magnetosphere has been approximated by an 
axisymmetric shape, it follows that the resulting field lines for this 
case are also confined entirely to the same plane. Results of two such 
calculations are shown in Figure 3* The magnetic field in the incident 
stream is inclined at 90° to the direction of the velocity vector in the 
left portion of Figure 3, and 45° in the right portion. The correspond- 
ing results for 0° inclination are not presented, but can be visualized 
easily because the field lines for that case are aligned everywhere with 
the streamlines, and the field strength ratio is proportional to 

pv/p V . 

These results clearly shew how the magnetic field lines bend dis- 
continuously as they pass through the bow wave at any angle except a 
right angle, and then curve in a continuous manner throughout the entire 


rer;ion between the bow wave and the maeneto sphere . The discontinuous 
bead at the shock wave is, moreover, always in the direction that pre- 
serves the sign of the tangential component of the field, as is required 
for all physically relevant hydromagnetic shock waves. It may he seen 
that the field lines illustrated in Figure 3 are all draped around the 
nose of the magnetosphere. Outside of this plane, however, the field 
lines drift past the nose with the flow and deform into three-dimensional 
curves. The strong constraint imposed on the magnetic field by the 
stagnation point is thus greatly reduced, and the field lines may be 
anticipated to remain much straighten than illustrated for the plane 
of symmetry. 

The choice of the value 2 for the ratio of specific heats employed 
in calculating Figures 2 and 3 is usually Justified by reference to the 
presumed two -degree -of -freedom, nature of the interactions of charged 
particles in a magnetic field. This argument weakens, however, when 
consideration is given to the irregular character of the magnetic fields 
observed in space, particularly downstream of the bow shock wave. In 
fact, the whole concept of applying hydromagnetic theory to the flow of 
solar plasma around the magnetosphere involving as it does the assumption 
of an isotropic pressure appears more consistent internally if the par- 
ticles are considered to behave as if they have three rather than two 
degrees of freedom and j is equated to 5/3- The effects on the veloc- 
ity and temperature contours of changing / from 2 to 5/3 are shown in 
Figure k. 

Although the values of 8 and 8 . 7 I employed in Figures 2, 3 , and 4 
for the free-stream Mach number are well centered in the range of values 


t 
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ti' be expeeted In 1 he incident solar wind, considerable variation is known 
to occi;r. In order to illvistrate the magnitude of the effects to be antic- 
ipated, results are presented in Figure 5 f’oi’ Mach numbers of 5 and 12 
for a gas with y - 5 / 3 - They show that the bow wave recedes from the 
magnetosphere as the Mach number diminishes and as 7 increases. The 
change is small in the ranges covered, however, as is the change in the 
velocity contours. Temperatures in the flow field depend strongly on 
both Mach number and /, however, with higher values associated with 
higher Mach number and larger 7. 

A useful q^uantity for characterizing the location of the bow shock 
wave is the standoff distance A at the nose of the magnetosphere. This 
distance has been shown in an aerodynamic context by Seiff (1962) and 
Inouye (1965) to be very nearly proportional to the density ratio across 
the nose of the bow wave for a wide range of values for 7 and The 

degree to which this simple empirical result is able to represent the 
standoff distance is illustrated in Figure 6. 

In conclusion, we have presented a sample of the type of results 
that can be obtained using the continuum equations of magnetohydrodynamics 
and gasdynamics. They augment the results for the location of the bow 
wave given by Spreiter and Jones (1963) which has been widely used in the 
interpretation of data from IMP -1 satellite and Mariner 2 and space 
probes. It has been reported in this syraposi-um by Goon that data from 
the Vela satellites are in at least qualitative agreement with these 
results. With respect to such comparisons, it should be observed that 
our calculations only yield information regarding the average bulk prop- 


erties of the flow, whereas the usual plasma probes and magnetometers 
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ol.u'erve a jmich more microscopic and instantaiaeous view of the plasma. 
With respect to the latter, it should also be recalled that substantial 
quantitative discrepancies still remain between the values for bulk 
parameters deduced from readings of plasma probes of different design 
even when flown on the same spacecraft. These differences are certain 
to he resolved as better understanding develops, and we present these 
results at this time in order that they be available for comparison 
with observational data as the latter becomes increasingly available. 
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FIGURE CAPTIONS 


Figure 1.- Streamlines and wave patterns for supersonic flow past the 
magnetosphere; = 8.7I, J ^ 2. 

Figure 2.- Contour maps of constant velocity^ mass flux, density, and 
temperature for supersonic flow past the magnetosphere; Moo 8-71^ 

7 - 2. 

Figure 3>“ Magnetic field lines in plane of free-stream velocity and 

magnetic field vectors for supersonic flow past the magnetosphere; 

8.71, 7 = 2. 

Figure 4.- Effect on velocity and temperature contours of changing ratio 
of specific heats from 2 to 5/S; Ho “ 

Figure 5*- Effect on velocity and temperature contours of changing Mach 
number to 5 and 12; 7 = 5/3* 

Figure 6.- Variation of standoff distance with density ratio across the 
nose of the how shock wave. 





Fig. 2. Contour maps of constant velocity^ mass flux, density, and tem- 
peraoure for supersonic flow past the magnetosphere: 1-^ = 5-7F; /' = 2. 


MAHONAl AEKC^UItCS AiUi-SPACl AOAUNlS^qUTlOH 
AMES 8I5CARCH CXNTR, MOFFETT FIELD. CALIFORNIA 


1 _ 



NJtTTOWAl A Et ONAU T tC S AW SP*a- A0WWlST1wn0»« 
AMES RESEARCH CEKTER. MOFFEH ESaD, CAL^ORNIA 



y=2 r-6/3 




ty and teiiperature con-cours cf changing ratio 


N*TTONf*l *F«ON*UTJC5 JWO-SFACt AOIWNTSUTATTON- 
AMES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 


M = 5 

-r 3 


M = 12 
T 3 



^02T. or: 




ar.d temperature 

5 ani 12; > 


con uo^ors 


-1 


changing Mach 


-KAnOKA-L AERONAUTKS ANa^SPAq ADMtNtSTOAIION 
AMES RESEARCH CEKTER. MOFFFn FIELD. CALIFORNIA 



(y-i)M|4 2 
/^b'(y+i)M|) 

DENSITY RATIO ACROSS NORMAL SHOCK WAVE 


maci 


iistaac 


NATKJNAL XaONAUTjCS AND Sf ACf AOMINISflA^ 
AMES RESEARCH CENTER. MOFFETT F»EID, CAUFORNIA 


